Phenomenological model of the Kaonic Nuclear 
Cluster K~pnn in the ground state 

A. N. Ivanov 1 ' 3 * 1 ", P. Kienle 2 ' 3 * J. Marton 3§ , E. Widmann 31 

February 9, 2008 



1 Atominstitut der Osterreichischen Universitdten, Technische Universitdt Wien, 
Wiedner Hauptstrasse 8-10, A-1040 Wien, Osterreich, 
2 Physik Department, Technische Universitdt Miinchen, D-85748 Garching, Germany, 
3 Stefan Meyer Institut fur subatomare Physik, Osterreichische Akademie der 
Wissenschaften, Boltzmanngasse 3, A-1090, Wien, Osterreich 

Abstract 

A phenomenological model, proposed for the Kaonic Nuclear Cluster (KNC) 
K~pp (or *>H), is applied to the theoretical analysis of the KNC K~pnn (or S°). 
The theoretical values of the binding energy and the width are equal to e^o = 

— 197MeV and T s o = 16MeV. They agree well with the experimental data e^ p = 

- 194.0i;i;! MeV and T 7? < 21MeV (PLB 597, 263 (2004)) and the theoretical 
values e$o = — 191 MeV and T^o = 13MeV, obtained by Akaishi and Yamazaki 
within the potential model approach. 
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1 Introduction 

Recently [1] Suzuki et al. have announced the experimental discovery of the bound state 
K~pnn, named S°(3115), with quantum numbers I{J P ) = 1(§ ), the binding energy 
4 X P = - 194.01^ MeV and the width constrained by Vg? < 21MeV. The possible 
existence of such a Kaonic Nuclear Cluster (KNC) as well as of the KNC K~pp (or 
which has been observed by the FINUDA Collaboration [2J, has been pointed out by 
Akaishi and Yamazaki within the potential model approach [3] [5]. According to [Tj, the 
main decay channels of the KNC 5' (3115) are S° — > H~pn, S° — > T,~d and S° — > TPnn, 
where d is the deuteron. The decay channel S° — > A°nn is suppressed. 

In this paper we extend the phenomenological quantum field theoretic model for the 
KNC ^H, proposed in [6], to the description of the KNC 5' (3115). The binding energy 
and the width of the KNC S° (3115) we define as [6] 



where dr is an element of the phase space of the K~pnn system, $50 is the wave function 
of the KNC S^^llS) and M(K~pnn — > K~pnn) is the amplitude of K~pnn scattering. 

In our model for the calculation of the scattering amplitude of the constituents of the 
KNC we use the chiral Lagrangian with a non-linear realization of chiral SU(3) x SU(3) 
symmetry and derivative meson-baryon couplings [7], accepted for the analysis of strong 
low-energy interactions in ChPT [8] (see also [9], [10]). We calculate these amplitudes in 
the tree-approximation to leading order in the large Nc and chiral expansion, where Nc 
is the number of quark colour degrees of freedom [TT]. According to Witten the large 
Nc expansion is equivalent to the heavy-baryon approximation, applied in ChPT [8]-[TU] 
to the analysis of low-energy meson-baryon interactions (see also [T2], [13] ) . 

In our model the main part of correlations between constituents of the KNC are 
described by the wave function of the KNC. 

For the construction of the wave function $ 5 o one has to specify the structure of 
the KNC S' (3115). The simplest structure of the bound K~pnn system with quantum 
numbers I{J P ) = 1(| + ) and the decay mode S° — > A°nn suppressed is (K~p) I=0 (nn)fz l 1 , 
where the K~p pair is in the state with isospin 1 = 0, i.e. (K~p)i =Q , and the nn pair is 
in the P-wave and spin-triplet state with isospin 1 = 1, i.e. (nn)^. 

The paper is organised as follows. In Section 2 we calculate the binding energy of the 
KNC S' (3115). We get 650 = — 197MeV. According to [6], the dominant contribution 
comes from the Weinberg-Tomozawa term for K~pnn scattering. In Section 3 we calculate 
the partial widths of the decay modes S° — > S~pn, S° T,~d, S° — > TPnn and S° — > 
A°nn and the total width of the KNC 5°(3115). We get r<?o = 16 MeV. We show that the 
dominant decay mode is S° — > Y,~pn with the partial width r(S' — > Yi~pn) = 15 MeV. 
We calculate the invariant-mass spectrum of the S° — > Y,~pn decay. In the Conclusion 
we discuss the obtained results. In the Appendix we define the wave function $50 of the 
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KNC S°(3115). 
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Figure 1: Feynman diagrams giving in the tree-approximation contributions to the am- 
plitude of the reaction K~pnn — > K~pnn to next-to-leading order in the large Nc and 
chiral expansion with respect to the Weinberg-Tomozawa term. 

2 Binding energy of the KNC 5° (3115) 

According to [6J, the binding energy of the KNC S' (3115) is given by 

- * ^ «) 1/2 = 197M eV, (2.1) 
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where F n = 92.4 MeV is the PCAC constant of pseudoscalar mesons. Other parameters 
co>i and Qi are related to the wave function $ 5 o of the KNC 5' (3115) and defined 
in the Appendix. The theoretical value agrees well with the experimental one e^ p = 
- 194.01 44 MeV. The binding energy is defined by the contribution of the Weinberg- 
Tomozawa term. This term gives the dominant contribution to the real part of the am- 
plitude of elastic K~pnn scattering to leading order in the large Nc and chiral expansion. 

The correction of next-to-leading order in the large Nc and chiral expansion to the 
Weinberg-Tomozawa term is defined by the diagrams in Fig. 1. The calculation of dia- 
grams in Fig. 1 is carried out using the chiral Lagrangian with a non-linear realization of 
chiral SU(3) x SU(3) symmetry and derivative meson-baryon couplings [6j. For Nc = 3 
such a correction makes up a few percent of the Weinberg-Tomozawa term. 



3 Width of the KNC 5° (3115) 

In this section we calculate the partial widths of the decay modes S° — > Y,~pn, S° — > H~d, 
S° -> S°nn and S° -> A°nn and the total width T s o of the KNC S° (3 11 5). For the total 
width we get the value r^o = 16 MeV, agreeing well with the experimental data [1] 
T s o < 21 MeV and the potential model prediction T s o = 13 MeV jl] (see also [5]). 

3.1 The S° — > YTpn decay mode 

The Feynman diagrams of the reaction K~pnn — > S~pn are depicted in Fig. 2. In the 
tree-approximation these diagrams give the main contribution to the amplitude of the 
S° — > H^pn decay to leading order in the large Nc and chiral expansion. The calculation 
of these diagrams is carried out using the chiral Lagrangian with chiral SU(3) x SU(3) 
symmetry and derivative meson-baryon couplings (see Appendix B of Ref.|[Sj). The partial 
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Figure 2: Feynman diagrams of the amplitude of the reaction K~pnn — > H~pn giving in 
the tree-approximation the main contribution to leading order in the large Nc and chiral 
expansion. 

width of the 5° — > E~pn decay is equal to 

■ :; ' 2<) •' '"' -/'.^.V 2 9 gl NN /2/i 3 



7T 



/ 327T 2 rriKm^F^ \1 



xe 2 f(e) = 15MeV, (3.1) 

where g^Aw = Qa^nIF-k = 13.3 is the 7riViV coupling constant [HJ [15], s = m s o — 
2m at — ms = 44 MeV is the energy excess calculated for m^ = 940 and ms = 1193 MeV, 
\i = 2mArms/(2mAr + ms) = 730 MeV and f(e) = 0.78 is the contribution of the phase 
volume of the T,~pn state. The parameters /ij, u>i for (i=l,2,3) and VL\ of the wave function 
$50 of the KNC S^^llS) are defined in the Appendix. 

3.2 Invariant— mass spectrum of the S° —> Y<~pn decay 

Following [16] we calculate the invariant-mass spectrum of the S° — > Yi~pn decay. For 
this aim we use the following variables [16] 



X = m 2 np = (E n + E p ) 2 - (k n + k p ) 2 , Y = m% = (E s - + E p ) 2 - (A* + k p ) 2 , (3.2) 

which have the meaning of the invariant-squared masses of the np and E~p pairs, respec- 
tively. At the rest frame of the KNC S' (3115) in terms of the variables (X, Y) the kinetic 
energies of the neutron and the proton are given by 

rp _ (mso ~ m N ) 2 - Y T _^ (mso -m s ) 2 + (m 5 o -m N ) 2 -X-Y 
2mso ' v 2mso 

The invariant-mass spectrum or the Dalitz density distribution of the S° — > E~pn decay 
is defined by 

3 d 2 T(S° -» Trpn) / filufQ^ 1 / 2 / fi 2 uJ 2 \^ 2 ^3UJ 3 \ 5/2 9 gl NN m^m s o 
ms ° dXdY \ 7T 3 J \ 7T J \ 7T J 128tt 2 m K m%Fl 

x [A{X, Y) + B{X, Y) + C(X, Y) + D(X, Y)] 2 = 
= 10[A(X,Y) + B(X,Y) + C(X,Y) + D(X,Y)] 2 - (3.4) 



In the (X, y)-plane the Dalitz domain is bounded by the curve with A min = 4m 
A max = (m s o - m E ) 2 , F min = (m s + m N ) 2 and F max = (m 5 o - m N ) 2 [16J. The invariant- 
mass spectrum of the S° — > H^pn is represented in Fig. 3. 
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Figure 3: The invariant-mass spectrum of the S° — > T decay in the variables X = m 2 
and Y = m| [16] . 
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3.3 The 5° -> Trd decay mode 

The partial width of the S° — > X _ <i decay is equal to [T7j 



r(s c 



x 



-> rrd) 



7T" 



-(A + 5 + C + D) 



2 np np 



lMeV, 



(3.5) 



where the amplitudes A, B, C and D are defined by the corresponding Feynman diagrams 
in Fig. 2, Qy,ci = 266 MeV is a relative momentum of the Trd pair, = y/fPiFy/ ( r ) 2 +k 2 ) 

is the wave function of the deuteron in the ground state in the momentum representation 
and 7 = 1/Rd = 46 MeV [18]. The amplitude of the np scattering is defined according to 
(see also [T7]) and taken in the effective range approximation 



k cot S np (k) 



1 1 t , 
i r l u 

™ ' np n i 



(3.6) 



np 



where S np (k) is the phase shift of the np scattering in the 3 Si state, 



np 



(5.424±0.004) fm 

and r^p = (1.759 ± 0.005) fm are the S-wave scattering length and effective range of np 
scattering in the 3 Si state [18J. 

The main contribution to the integral comes from the relative momenta k ~ 100 MeV, 
and 94 % of the value of the integral are concentrated in the region < k < 200 MeV. For 
the relative momenta k = 100 MeV and k = 200 MeV the phase shift of the np scattering, 
defined by the effective range approximation (13.61) . is equal to 5 np = 85.3° and 5 np = 54.6°, 
respectively. These values agree well with the SAID analysis of the experimental data on 
the np scattering in the 3 Si state [20]: 5 np = 85.84° and 5 np = 48.84°, respectively. 
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Figure 4: Feynman diagrams of the amplitudes of the reactions K'pnn hPnn and 
K~pnn — > E°nn giving in the tree-approximation the main contribution to leading order 
in the large Nc and chiral expansion. 

3.4 The S° — ► E°nn and 5° — ► A°nn decay modes 

The amplitudes of the S° — > and S° — > A°nn decays are defined by the integrals 

of the amplitudes of the reactions K~pnn — > TPnn and K~pnn — > A°nra with the wave 
function $50 of the KNC 5' (3115). Feynman diagrams of the amplitudes of the reactions 
K~pnn — > S°nn and K~pnn — > A°nn are depicted in Fig.4. In the tree-approximation 
these diagrams give the main contribution to leading order in the large iV^ and chiral 
expansion. The contributions of the amplitudes of the reactions K~pnn — > TPnn and 
K~pnn — > A°nn, calculated to leading order in the large Nc and chiral expansion, to the 
amplitudes of the S° — > S°nn and S* — > A°nn decays vanish. This means that the decay 
modes S* — > TPnn and S* — > hPnn are suppressed to leading order in the large Nc and 
chiral expansion. 

3.5 Width of the KNC 5° (31 15) 

The width of the KNC S' (3115) is defined by the sum of the partial widths of all modes. 
Since the dominant modes are S° — ► T~pn and S° — > T~d, the width of the S°(3115) is 
equal to 

r 5 o = 16MeV. (3.7) 

This agrees well with the experimental constraint r^, p < 21 MeV and the theoretical value 
T^o = 13 MeV, predicted by Akaishi and Yamazaki within the potential model approach 
[1] (see also [5]). 

4 Conclusion 

The phenomenological quantum field theoretic model of the KNC K~pp (or g-H), pro- 
posed in [6], is extended to the description of the KNC S' (3115) with quantum numbers 
I(J P ) = 1(| ) p]. The wave function of the KNC S^SHS) is taken in the oscillator 
form. The frequencies of oscillations are related to the frequency u, describing longitu- 
dinal oscillations of the 7^~-meson with respect to the pp pair in the KNC [6]. The 
value of the frequency uj = 64 MeV has been calculated in [6] by fitting the width of the 
A°(1405) hyperon, treated as a bound K~p state ^H. 

We would like to make clear the title of our model namely "The phenomenological 
quantum field theoretic model" . On the one hand our model is "phenomenological" , since 
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we construct the wave function of the KNC, the parameters of which are fixed in terms of 
the experimental data on the width and mass of the A (1405) hyperon. On the other hand 
our model is "quantum field theoretic" , since to the calculation of scattering amplitudes 
of the constituents of the KNC we apply the chiral Lagrangian with chiral SU(3) x SU(3) 
symmetry and derivative meson-baryon couplings, ChPT and the large Nc expansion. 

The binding energy e^o = — 197.0 MeV and the width V s o = 16 MeV, calculated in our 
model, agree well with the experimental data e^, p = - 194.0+ \\ MeV and < 21 MeV. 
We have found that the decay mode S° — > TTpn is dominant, whereas the decay modes 
S° — > TPnn and S° — ► A°nn are suppressed. The calculation of the amplitudes of elastic 
and inelastic K~pnn scattering we have carried out using the chiral Lagrangian with the 
non-linear realization of chiral 577(3) x SU(3) symmetry and derivative meson-baryon 
couplings (see Appendix B of Ref. [6]). We have kept the contributions of leading order 
in the large Nq and chiral expansion [6] . In more detail the experimental analysis of the 
theoretical amplitude of the S° — > H~pn decay can be obtained through the measurements 
of the invariant-mass spectrum [16]. The theoretical invariant-mass spectrum of the 
S° — > £~pn decay is given by Eq. fl3.4p and represented in Fig. 3. 

We would like to accentuate that our theoretical predictions for the binding energy 
e s0 = - 197.0 MeV and the width I> = 16 MeV of the KNC 5° (3115) agree well with 
results, obtained by Akaishi and Yamazaki within the potential model approach [4] (see 
also [5]): e 5 o = - 191.0 MeV and T(S° -> E~pn) = 13 MeV. 

However one can show that in our model the nuclear matter density of the KNC 
S°(311B) is smaller by factor 3 than nso(0) = 1.56 fm -3 , predicted in the potential model 
approach [5]. Indeed, using the wave function of the KNC S°(3115) one can calculate the 
nuclear matter density at the center of mass of the K~pnn system. This gives 



The obtained result is larger by a factor 4 than the normal nuclear matter density no = 
0.14fm -3 , but it is smaller by factor 3 than n s o(0) = 1.56 fm~ 3 , predicted by Akaishi and 
Yamazaki. 

Our approach for the analysis of the Kaonic Nuclear Clusters K~p (or A°(1405)), K~pp 
and K~pnn can be applied to the description of the Kaonic Nuclear Clusters K~ppn with 
isospin 1 = 0, K~ppn and K~ppp with isospin 1 = 1, K~ NNNN and so on [5]. 



Appendix A: Wave function of the KNC 5°(3115) 



In the momentum representation the wave function $50 of the KNC S' (3115) is defined 




(4.8) 





7 



<5>nn{Q) 



( 



S 7T ^3/4 
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Q 2 \ 



2fi 3 u 3 



(A-2) 



where the nn pair is in the P-wave state, Yi m ($A, ^a) is a spherical harmonic, m = 0, ±1 
is a magnetic quantum number. The frequencies and the reduced masses are 



fix 

0J2 
^3 



3 /i 7KA/r Tr Sm^m^ 

c<j = 76 MeV , /ii = 

2 /i! ' mjf + 3m N 



420 MeV 



V^wi = 132 MeV, 




u 2 
— c<j = 51 MeV , \i2 = - m N = 
to 3 

- [JLuj = 30 MeV , u 3 = -m N 
2 V to 2 



630 MeV, 
= 470 MeV, 



(A-3) 



where cj = 64 MeV is the frequency of the longitudinal oscillation of the i^~-meson 
relative to the pp pair in the KNC |.H [6] and \i = 2niKmN / {mK + ^m^) = 391 MeV is the 
reduced mass of the K~(pp) system, calculated for tuk = 494 MeV and = 940 MeV. 
The value us = 64 MeV is obtained from the fit of the width of the A°(1405) hyperon, 
treated SIS db bound K p state [5J. 

The wave functions $ p ( nn )(fc) and <& nn {Q) are normalised by 



d 3 q 



1$ 



p(nn) 



This gives 
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(A-4) 



(A-5) 



The wave functions ^k(q), & P (nn)(k) and § n n(Q) describe oscillations of the i^~-meson 
relative to the pnn system, the proton relative to the nn pair and the neutrons in the nn 
pair, respectively. The K~pnn system with the wave function (1 A- 1 [) possesses positive 
parity. 

Since the total angular momentum of the KNC S' (3115) is J = § [TJ HI [5] , the wave 
functions ^jm{Q) of the angular momentum of the KNC S' (3115) are defined by 
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where XjTli with a = l,2,p are the spinorial wave functions of the neutrons and the 
proton, respectively. 
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